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The problem of an electromagnetic wave scattering by a slab with two rough bound- 
aries is solved by a small-perturbation method under the Rayleigh hypothesis. In 
order to obtain a perturbative development, we use a systematic procedure which 
involves integral equations called the reduced Rayleigh equations. Then we will show 
for a dielectric slab deposited on a silver film that the backscattering enhancement 
can be produced by guided waves which interact with the two rough surfaces. 
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1 I. INTRODUCTION 

Scattering of electromagnetic waves from multilayer structures is a phenomenon which is of interest in 
^ . many area of physics such as remote sensing or optical industry, where for example metallic surfaces have 
a dielectric coating. An extended review of experimental andj-theoretical woritSpjOn optical multilayers 
can be found iru, and references on recent papers are given ina. Many worksalaffl deal with the small 
perturbation method to investigate the behaviour of multilayer structure. However, due to calculations 
complexity, the analytical results are given only to first order in the rms- height of each rough boundary, 
thus the interaction between the rough surfaces can only be taken into account if the rough surfaces are 
correlatedu in the mean procedure. Our purpose is to show how light can interact with several rough 
surfaces, to this end, we have choosen the most simple system depicted in Fig. [j], where we have three 
regions with different permitivities separated by two rough surfaces. The calculations [have been made 
under the hypothesis of the small perturbation method, initially developped by RiceQ'Q. Due to the 
great complexity to derive high orders of the perturbative development, Rice's original method is rather 
difficult to apply. One way to overcome these difficulties isJjO use integral equations called the reduced 
Rayleigh equations, they were first obtained by Celli et aE, and later have been generalized in Ref.Q 
£j | to take into account upward fields in the slab medium. When combining these equations we obtain an 
• • ■ integral equation where only the scattering matrix of the whole structure has to be determined, next 
we have developed a systematic method to calculate the perturbative development. We have shown by 
means of numerical simulations, how different mechanisms, responsible of the enhanced backscattering 
occur. This phenomenon which has been predictecOOE3t3 and observedtil in the case of a random 
metal surface, manifests itself as a well-defined peak in the retroreflection direction through the angular 
dependence of the intensity related to the diffuse component of the scattered field. For a metallic surface, 
the phenomenon is produced by the interference of waves which excite surface plasmon polariton along 
a certain path and then follow the spne path but in the reverse direction. For a dielectric bounded by 
metallic plate with one rough surfacec3, the enhanced backscattering is produced by a similar mechanism 
where the surface plasmon polariton is replaced by a slab guided wave. In this case the incident wave 
excites in a first time a guided mode due to surface roughness, and then the roughness transforms 
the surface Wjuie into a bulk wave. Furthermore, if the slab supports several guided modes, recent 
investigationsoiij have shown the presence of additionnal peaks, called satellite peaks, in the angular 
distribution of the incoherent intensity. When the slab has two rough boundaries, the enhancement of 
backscattering can be produced by two kinds of interaction. The classic one, where the wave is scattered 
two times by the same boundary, and a new one (see Fig. where the wave is transmitted at a point A 
of the surface without being diffused, then the wave is scattered by the second rough surface at a point 
B and by the first one in C . If the wave follows the same path in the reverse direction and excites the 
same guided mode, then phase difference between the two path is Acj) = tbc ' (fc + where tbc is 
the distance between B and C. Thus, in the anti-specular direction (k = —ki), this phase difference is 
independent of the random position of B and C which produce the backscattering peak. 

The outline of the paper is as follows. In Section [n] we present the system under study and a description 
of the rough surfaces statistics. In Section |TJ we define the plane wave representation of the electric 



field in the polarization basis and also the scattering matrix. Section IV is devoted to the calculation of 
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the incoherent cross-section as a function of the perturbative development. In Section [v| we derive an 
integral equation from the reduced Rayleigh equations, which provides a systematic way to obtain the 
perturbative development. The resulting expressions are numerically computed in the case of a dielectric 
slab deposited on a rough silver surface in Section [v|. The calculated incoherent intensity will show 
a narrow enhanced backscattering. Conclusions drawn from our results are presented and discussed in 



Section VII. In an Appendix, we collect matrices used in the calculations. 



II. THE RANDOM SURFACE 

The system which is considered in this paper is depicted in Fig. |[ 

The three regions are characterized by an isotropic, homogeous dielectric constant eo, ei and £2 respec- 
tively. The two boundaries are located at z = h\{x) and z = —H + h,2(x), x = (x,y), and these three 
media are separated by rough surfaces described statistically. In fact, we assume that h\{x) and h,2(x) 
are stationary, isotropic uncorrelated Gaussian random processes defined by their moments : 

<hi{x)> = 0, (1) 
<h i (x)h i {x') > = W i (x-x / ), (2) 
< hi(x) h 2 (x') > = (3) 

where i = 1,2, and the angle brackets denote an average over the ensemble of realizations of the function 
h\(x) and hiix). In this work we will use a Gaussian form for the surface-height correlation function 
Wi(x) and W 2 (x): 

W i (x)=a* exp(-x 2 /l*), (4) 

where cti is the rms height of the surface hi(x), and li is the transverse correlation length. In momentum 
space we have: 

<hi(p)> = 0, (5) 

<h l {p)h l (p l ) > = (2ir) 2 5(p+p')W l (p), (6) 

<h 1 ( P )h 2 (p')> = 0. (7) 



with 



Wi(p) = J d 2 xW l (x) exp(-ip-x) (8) 
= na 2 l 2 exp(-p 2 / 4 2 /4). (9) 



III. THE SCATTERING MATRIX 

We suppose that the slab is illuminated from the media by an electromagnetic wave of pulsation to. 
In the following we will omit the time dependence exp(— iu>t). The field E° in the media can be written 
as a superposition of an incident and scattered fields : 

j^-^E s (p) exp(ip ■ x + ia Q (p) z) . (10) 



where 



a Q (p) = (e K 2 ~p 2 )i , (11) 
E\p Q ) - £MPo) e° v -{p ) + E l H (p Q ) e H (p Q ) (12) 
E s (p) = E° v {p) e° v + (p) + E s H {p) e H (p) . (13) 
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The subscript H refers to the horizontal polarization (TE) and V to the vertical polarization (TM), and 
are defined by the two vectors: 



e H (p) = e z x p, (14) 
x0±,_^ _ j_ < *o(p) ~ _ [jgj 



4 ± (p)=±^P-^e z , (15) 



where the minus sign refers to incident wave and the plus sign to the scattered wave. It has to be noticed 
that the vector E s (p) and E l (p ) are expressed in different basis due the fact that e.y^(p) and e]^{p) 
depend on p. In medium 1, we have a similar expression namely: 

/c]l^P f d^p 

-^-^E 1 -{p)exp(ip-x-iai(p)z) + J -^-^E 1+ (p) exp(ip • x + i£*i(p» , (16) 

where 

a 1 (p)={e 1 K^-p 2 )i. (17) 
The field E 1 ~ is decomposed in the basis (e.y(p), ejj{p)), and E 1+ in the basis (e.y(p), eu(p)) with 

e H {p) = e z xp 1 (18) 

(19) 

We now introduce the definition the scattering matrix : 

E s (p)^R(p\p )-E\p ), (20) 

where R(p\p ) is a two dimensional matrix which can be written in the following form: 

pr^lr, \ - f R w(p\Po) Rvh{p\p )\ 
R[P\Po) ~ [r hv ( p \ Po ) R HH (p\p )J ■ 

IV. DIFFUSE CROSS-SECTION 

In a previous worlH we have defined a new product between two dimensional matrices in the form : 
fvv fvH \ ~, ( 9vv 9vh 



f®g= 7 Y © r yv " (2i) 

V/ffV JHHj \9HV 9hhJ 

( fw9vv fvH9v H Re(fvvgvn) -Im(fvvg VH ) 

fHvg* HV fHHg*HH Re(f HV g* HH ) -Im(f H vg V H) 

2Re{f vv g* HV ) 2Re{f VH g* HH ) Re{fvvg'v V + fHvg VH ) -Im(f vv g* HH - fvHg* HV ) 

\ 2Im{fw9*Hv) 2Im (fvHg* H H) Im (fvvg vv + fHvgvn) Re(fvvg* HH - fvHg* HV ) 

which allows to write the incoherent Muller matrbS in a condensed expression: 



M nc ° h (p\p ) = ° 2 [< R(p\ Po ) R(p\ Po ) >-< R(p\p ) > < R(P\P ) >] ■ (22) 



incoh f _\_ s KqCOS 2 8 f — f 

Wf 

Furthermore, we define the generalization of the classical bistatic coefficientO also called mean differential 
coefficientlia by : 



~f" lc ° h (p\p ) = -^—M incoh ( p \p ) (23) 
A cos 0q 



where A is the area of the illuminated surface and 9q the incident angle (see Fig. ^). In this paper we are 
interested by the perturbative development of the scattered fields as a function of the surface elevations 
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h\ and h 2 - In a perturbative expansion of the scattering matrix, the terms which contain an expression 
like h™ will be denoted R^ nm \ so that the perturbative development of -R becomes: 



,(00) 



j(10) 



5(01) 



j(20) 



j(21) 



l(12) 



l(22) 



j(30) 



j(03) 



R = R y ""' + R y "' + R" + R" + R"' + R"' + R~' + R"' + R y ""' + IV""' + ... (24) 
With this decomposition and using Eqs. (0,0) with the fact that hi and h 2 are Gaussian random processes 



< h 2p+1 (x) > = , pa positive integer 

< h 2 2 p+1 {x) > = 



(25) 
(26) 



the incoherent bistatic matrix will be given by the contribution of three terms: 



i mcoh (p\Po) = r: coh {p\Po) +ir on (p\Po) +iz on (p\p ) , 



—incoh 1 



—incoh 1 



(27) 



where 



A(2tt) 2 cos6» 



< R iW) R m > + < R m R m > + < R m R iW) > 



(28) 



corresponds to the incoherent bistatic matrix for the slab where only the upper surface has a roughness 
(fi2(x) = 0), and its expansion is made up to order four in the rms-height elevation ax- Similarly, 



—incoh 
Id 



(p\Po) 



Kl cos 2 e 
A(2tt) 2 cost 



<R m QR {01) > 



<R m QR m > + <R i03) QR m > 



(29) 



is associated to a system where only the bottom surface is rough (hi(x) = 0), where also the perturbative 
development is made up to order four in a%. The last term ^^d contains terms which describe the 
scattering process between the two rough surfaces where only the leading terms are retained: 



lud COh (P\Po) 



Kl cos 2 e 
A{2it) 2 cos 6»o 



-=(10) -=r(12) 

< R K ' QR > 



—(01) t=(21) 

< R y ' R y ' > - 

—(11) -=(n) 

< ,R V ' JT ' > 



-=(12) -=(10) 

<R y ' R y ' > 



<R {21) QR m > 



(30) 



More precisely all terms up of order a\ erf are included. If the value of <j\ and o 2 are of the same order 
of magnitude, these terms will be comparable to the order four in the expressions (|28[|29] ) ■ Thus we can 
suppose that the following terms in the expansion d30|), which are of order u\ erf, o\ erf, o\ cr| will be 
negligible compared to those retained in Eq. (^0|). However, due to the complexity of the perturbative 
development, these terms of the sixth order have not been calculated. In the following section, we will 
show how the perturbative development can be put in the following form: 



-R (10) (p|Po) = "o(Po) ^ (10) (pIPo) Mp - Po) : 
-R (01) (p|Po) = «o(p ) ^ (01) (p|Po) MP - Po) ■ 



(31) 
(32) 



R {11) (p\po) = MPo) J J^TZ 



-^ (11)12 (p|Pi|Po)MP - Pi) Mp-Po) 
+X (11)21 (p|p 1 |p )/i 2 (p -Pi) hxip-po) , 



(33) 



# (20) (p|Po) = "o(Po) J j^k X m '(pipxipo) hx(p - fti(pi - p ) , 
^R ( ° 2) (p|Po) = "o(Po) J j^X m {p\p l \p )h 2 {p-p l )h 2 {p l -p ), 



(34) 
(35) 



-R (21) (p|Po) = «o(Po) 



d 2 Pi d 2 p 2 
(2tt) 2 (2tt) 2 



' (21)112 

x (p\Pi\p 2 \Po) h i(P ~ Pi) h i(Pi ~P 2 )h 2 (p 2 -Pa) 

(21)121 

+ X (p|Pl|P 2 |P ) h l(P -Pl) h 2(Pl -P 2 ) ^l(P 2 -Po) 
(21)211 

+X " (p|Pi|p 2 |Po)^2(p -Pi)fti(Pi -P 2 ) hi(p 2 -Po) 



(36) 



i2 (12) (p|p )^a (Po) 



d 2 Pi d 2 p 2 

(2tt) 2 (2tt) 2 



(12)221 

X (p|Pl|P 2 |P ) ^(P -Pl)/l 2 (Pl -p 2 ) ft l(P2 "Po) 
(12)212 

" (p|Pi|p 2 |Po) /l 2(p-Pi)/ii(Pi -p 2 ) h 2 {p 2 -p ) 

(12) 122 

\-X " (p\Pi\p 2 \Po) h i(P -Pi)h 2 (Pi -P 2 ) ^ 2 (P 2 -Po) 



(37) 



-R (30) (p|Po) = ao(Po) 
-R (03) (p|Po) = "o(Po) 



d 2 Pi d 2 p 2 —(30) 



(2tt) 2 (2tt) 2 
d 2 Pi d 2 p 2 —(03) 



X (p|Pi|p 2 |Po)^i(p -P\)hi{Px -p 2 )MP2 -Po) j (38) 



(2tt) 2 (2tt) 2 



-X (p|Pi|P 2 |Po)MP -p 1 )h 2 (p 1 -p 2 )h 2 {p 2 -Po) ■ (39) 



In these expressions we have added superscripts in some terms to indicate the order of apparition of 

(21)1 2 1 

the functions h\ and h 2 . For example, in X the superscript 121 indicates that it is the coefficient 

associated with the product hi(p—p 1 ) h 2 (p 1 —p 2 ) hi(p 2 —p ). When we combine Eq. ml) , 5(0) = A/(2tt^ 2 
and the previous development, we obtain the following expression for the quantities (^8[]3 



r: coh (p\p ) = 

iT coh (p\p ) = 

rz coh (p\p ) = 



A o 


cos 2 9 cos 


% 




(27T) 2 






cos 2 9 cos 


% 




(27T) 2 






cos 2 9 cos 


% 


(27T) 2 



7 (10 - 10) (p|po) + 7 (20 - 20) (p|po)+7 (30 - 10) ' 



WOl-Ol). . . T (0 2 -02) -=(03-01), 
I (p\ Po ) + I (p\ Po ) + I 



-=(12-10) . . -=(11-11), , , -=(21-01), 

7 (p|Po)+J (p|p ) + 7 



mm 




(p\Po) 


(40) 


(p\Po) 


(41) 


(p\Po) 


(42) 



where 



/ (1O_1O) (p|Po)-^ 1 (p-po)X (1O) (p|po)0X (1O) (p|p o ) 



-=(20-20) 

I (p\Po) = 







d 2 Pl 
(2^) 2 



Wi(p -Pi)Wi(p x -po)X (20) (p|p 1 |p ) 



X (p|p 1 |p )+X (p|p+p -PilPo) 



I {3 °- 10 \p\Po) =Wx(p- Po) k (10) (p|p ) X (30) (p|p ) + X (30) (p|p ) X (10) (p|p ) 



(43) 

(44) 
(45) 



J (01 ~ 01) (p|Po) = W 2 (p - po) X (01) (p|p ) X (01) (p|p ) 



-=(02-02). 

I iP\Po) 



d 2 Pl 



>^( 02 ), 







2 w 2 (p- Pl )w 2 (p 1 - Po )x (pIpiIpo) 

-vf(° 2 )/ i i \ I -^(02). 



(271)2 

X v ""'(p|Pi|Po) 



(p|p+Po -PilPo) 



I ( ° 3 ~ 01) (p|Po) = Wa(p - Po) k (01) (p|Po) © ^ ( ° 3) (p|Po) + ^ ( ° 3) (p|Po) © * (01) (P|P ) , 



(46) 

(47) 
(48) 



() 



I (12 10) (p|Po) = ^i(p-Po) 



X (12) ( P |Po) X il0 \p\ Po ) + X {W \p\p Q ) X (12 \p\p ) 



i {21 ~ 01) (p\po) = w 2 ( P Po ) \x (21 \p\ Po ) x (01 W ) + X m (p\p ) X (21) (p\p ) 



i OpIpo) 



d 2 Pl 

(2tt)2 



Wi( P - Pi) W2(Pi - Po) ^ (U)12 (pIPiIPo) 



Q(x (n)12 (p\p 1 \p ) + X (im (p\p+p -p 1 \p ) 
+ W 2 (p-p 1 )W 1 ( Pl -p )X {11]21 (p\p 1 \p ) 



© (x (n)21 ( P | Pl | Po ) 



X (11)12 (p\p + p -PilPo) 



(49) 
(50) 



(51) 



with 



X (30) (p|p ) 



X (03) (p|p ) 



d 2 Pl 

(27T) 2 

d 2 Pl 

(2tt)2 



^ 1 (p 1 -p )X (30) (p|p |p 1 |p ) 

+^(p-p 1 )(X (30) (p|p 1 |p -p+p 1 |p ) + X (30) (p|p 1 |p|p ))] , (52) 
W 2 (p 1 -p )X (03) (p|p |p 1 |p ) 

+^2(p-p 1 )(X (03) (p|p 1 |p -p+p 1 |p ) + X (03) (p|p 1 |p|p ))l , (53) 



x (12 \p\ Po ) = 



X {21 \p\Po) 



d 2 Pl 

(2tt)2 

d 2 Pl 

(2tt)2 



W 2 (p - Pl ) (x (12)221 (p|Pilp|Po) + X (12)212 (p| Pl | Po - p +Pi|p ) 



fW 2 ( Pl -p )X (12)122 (p|p |p 1 |p ) 

W 1 (p~p 1 )(x {21)112 (p\p 1 \p\Po)^ 



x (21)121 ( P | Pl | Po - P + Pl | Po ) 



+W 1 {p 1 - p ) X (21)211 (p|p |p 1 |p ) , 



(54) 



(55) 



V. PERTURB ATIVE DEVELOPMENT AND REDUCED RAYLEIGH EQUATIONS 



In order to obtain the development (24) for the scattering matrix, a practical method is provided by 
the use of the reduced Rayleigh equations first obtained by Brown et alu for a single rough surface, and 
then extended inEI for a more general system. These equations are exact, under the Rayleigh hypothesis, 
and their main advantages are that one of the electric fields E l , E s , E 1 ~, E 1+ (see Fig. ||) of the problem 
has been eliminated. These equations are derived by taking linear combinations of the electromagnetic 
boundaries conditions at the first surface hi , where the Fourier transform of the fields has been introduced. 
In particular, we obtain the two following equations (see Eqs. (99-100) ina): 



{V ' P M l h + - a+ {u\p) ■ R(p\ Po ) ■ E\p ) + M* + '° (u|po) • E\p ) = 



(2vr) 2 

d 2 p 



(2^) 2 
where 



mI' 0+ (u\ p ) • R(p\p ) ■ E\p Q ) + M l h (u\p ) ■ E\p ) 



^yili^M ^+H, (56) 

(ei - eo) 
2(e ei)^ ai(ti) 



(ei - eo) 



M h (u\p) 



_ I{ba\{u) - aa (p)\u - P ) _ ib,o a 

ba\{u) — aao(p) U ^ 



E'-iu), (57) 



(58) 



with 



-r— -lb,0a . . , 

M (u\p) 



|tt|||| P || + abai(u) ao(p) u p — be§ Kq a\{u) (u x p) z 
ae? K a (p) (u x p) z (e ei)? K% it p 



(59) 
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and 



I(«|p) S /d>*e^Hp-*-i«M*)), 



(60) 



the symbols a = ± and 6 = ± in Eqs. (|39|-|60|) represent a given choice linked to the field propagation. 
In order to obtain a single equation for R(p\p ), we have to find a relation between E x ~ and E 1+ . To 
this end we already know an expression of the scattering matrix for a single rough surface separating two 
homogenous media of permittivity t\ and e 2 , which is translated along the z-axis to the height z = —H , 



and illuminated by a plane wave E 1 (p). This scattering matrix denoted £2 is given by 



R sei ,e 2 (p\Po) =exp(i(a 1 (p)+a 1 (p ))H)R sei , e2 {p\p ), (61) 

where R sei ^ 2 {p\p ) can be found in Refs.tMi3£jl The phase term in Eq. ( |6l| ) comes from the translation 
z = —H of the rough surface h 2 {x) (see RefsJUu). Thus we have the following relation: 



(62) 



Now combining Eq. ( p2| ) with Eqs. ( p6|f37| ), we obtain an integral equation where R(p\p ) is the only 
unknown: 



/ 



d 2 p 

(2tt) 2 



M h (u\p) + 



d 2 «i ai(u) -y-H 



M\ +fi -{u\p ) + 



(2tt) 2 ai(ui) 



«(p|Po) 



l-,0- 



(2tt) 2 ai(wi) 



R sei , e2 ( M l M i) ' M /> ' HPo 



Expanding /(a|p) in Eq. (|60|) in power of /ii: 



(63) 



/(a|p) = (2^) 2 5(p) - ia fc« (p) - ^ /if } (p) - ^- ft® (p) + ■ 



h[ n) (p) = J d 2 x exp(-zp • x) h%(x) , 



(64) 
(65) 



and using the perturbative development^ of R s e2 in power of h 2 



R s Ell£a (PlPo) = ( 27r ) ^(P - Po) X s ei , ea (P ) + «0(P ) * s ei,e a (PlPo) MP - P ) 



-«o(Po) 
-«o(Po) 



^fef.^CPlPllPo) h 2{P - Pl)h 2 (Pl - Po) 
d 2 p 1 d 2 p 2 -=H(3) / I I I \ 7 / \ , / \ , / \ 

(2T) 2 (2T) 2 X "i^(PlPilP2lPo) Mp - Pi)^2(Pi - P 2 )h 2 (p 2 - Po) , 



(66) 

we finally obtain the expansions (j2j^^-^3). The expression for the scattering matrix when only one 
rough surface is involved was given in RefH, here a small change have been made in the notations: 



R in °\p\Po) = R { u\p\Pa) , 

R ( " n \p\Po) = R { d\p\Po) 



(67) 
(68) 
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n being an integer ranging from to 3. For the others coefficients we have : 



xf(H)21 



^f(H)12 



(p\ Pl \ Po )=T (p)-U"(p)-X s : i ; e2 (p\p 1 )- -e-DM-X" '( Pl \p )+iS'( Pl \p 



x?(l°), 



(p|Pi|Po) =i-P + (p|Pi) • X y ^'( Pl \p ) , 



-=^(21)112 



X ' (pIPxIPalPbMP (plPx)-X 

l r 



vr(U)12. 



-^(21)121 
—(21)211 



(P1IP2IP0) 
ai(p)P (p|p 2 ) - a a {p 2 ) P + {p\p 2 ) 

(11)21, 



•^ (01) (P2|P ) 



(PIP1IP2IP0) =1^ (PlPi) ' ^ (P1IP2IP0) , 
(p|Pi|p 2 |Po) =T W (p)-U i0) (p)-X^l j (p\p 1 ) - [-e-D w { Pl 

, * ( £ i ^ £ o) , . yKio), , > 

:M ( Pl \p 2 )-X (p 2 \p ) 



^ (P1IP2IP0) 



2(e ei) 1 /2 

~ (aiCpJ^CPilPo) + ao(Po)S (PilPo) 



—(12)122 
-=^(12)212 



(P |Pi |P 2 |Po ) =i -P (p\Pi)- x (Pi |P 2 IPo 



-io 



r(°) 



-(11)12, 



(p|piIp 2 |po) =r (p) • t/ v (p) • ^ sei ,; 2 (p|Pi) • [-e- £>io(pi) ■ (Pi|P 2 |Po 

+ — — M ( Pl |p 2 )-X (p 2 |p ) 



2(e ei) 1 / 2 " 



(12)221 



(p|Pi|p 2 |Po) =^"(P) • C/ (0) (P) • [-X^CplpO-e - Oio(Pi) • ^"'"(PilPalPo) 
+ Xf e ( i 2) £2 (p|p 1 |p 2 ) • (-e - (p 2 ) ■ X m (p 2 \ Po ) +iS + (p 2 \ Po ) 



;H(1 



(11)21, 



where the matrices T 



7=(10) — (0) — - — ± -r— (10) W20) W01) ^r(02) -r— ff(l) ^H(2) 



(69) 
(70) 

(71) 
(72) 



(73) 
(74) 

(75) 

(76) 
(77) 



U 



P , D 



, X ea , X seue2 are respec- 



tively defined by Eqs.(71), (87), (114,115), (57), (104,105), (96,97), (60,61) in Ref.B, and e, S are given in 
Appendix |A|. 



VI. NUMERICAL RESULTS 



As an application of the previous formalism we consider a system made of a air-dielectric film whose 
dielectric constant is t\ = 2.6896 + i 0.0075(eo = 1), deposited on silver surface with e 2 = —18.3 + i 0.55. 
The vacuum-dielectric interface is a two-dimensionnal rough surface, whose parameters are o\ = 15 nm 
and = 100 nm. The dielectric-silver boundary is also rough and defined by a 2 = 5 nm and l 2 = 100 nm. 
The incident wave has an arbitrary polarization and his wavelength is A = 632.8 nm. With this set of 
parameters, the conditions of validity of the small-perturbation theory are satisfiedo namely: 



2tt 



1/2 



2tt 



1/2 



(78) 
(79) 



The thickness of the film is H — 500 nm and support two-guided wave polaritons for the (TE) polariza- 
tions at p\ E = 1.5534 K , and 

Pte = 1-2727 Kq, and three guided-modes for the (TM) polarizations at 
Ptm = 1-7752 Ko, Ptm — 1-4577 Kq and Pxm ~ 1-034X0- We have computed the incoherent bistatic co- 
efficient 7 mco,l (p|Po), where the integrals involved in Eqs. (44 47| , |5l| , p2[ - |55| ) are evaluated using Legendre 
quadrature. The results are shown in Figs. (§-0), where the incoherent bistatic coefficient is drawn as of 
function of the scattering angle 8 for two different angles of incidence and the incident wave is linearly 
polarized. In Fig. ^, the wave is normally incident and the scattered field is observed in the incident plane 
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(4> = 0°). The single scattering contribution on each surface, associated with the terms 1 °' > + I^ 01 01 \ 
is plotted as a dotted line, the double-scattering contribution I 20 20 + I 02 ° 2 ' ) + I 11 11 as a dashed 

-=(30-10) =(03-01) -=(12-10) =(21-01) . , , ... . . , 

line, the other terms I + I + 1 + I as a dash-dotted line, and the total contri- 

bution -f mcoh by the solid curve. We observed an enhancement of the backscattering which corresponds 
to the physical process in which the incident light excites a guided-mode through the roughness of the 
slab and then is scattered into a volume wave which is also due to the roughness effect. During the same 
process, the light can follow this path in the opposite direction where one possible configuration is shown 
in Fig. [|. These two paths can interfere constructively near the backscattering direction to produce a 
peak. As these paths are identical for the two waves under consideration, they should have the same 

degree of interaction with the rough surface, thus a term like J** 30 cannot produce the peak because 
the first wave interact three times with the upper rough surface while the second wave only one time. 

mi rr i r i —(20-20) =(02-02) —(11-11) ... . , .... , 

lhe effect comes only irom the terms I + I + 1 , which contain the paths indicated 

in Fig. |H However, it has to be noticed that these terms contain also paths that do not produce en- 

(20—20) 

hanced backscattering, for example in I we have the scattering process where the incident wave 

is only scattered one time by the upper rough surface but the scattering process is of order two in h\. 

ml , . „ , =(20-20) =(02-02) =(11-11) „ 

1 his is the reason why the terms I + I + I are not zero away from the anti-specular 

direction. In order to separate the different contributions to the backscattering peak, we have drawn the 

.=(20-20) =(02-02) =(11-11) , _ . _,. i-i , , , ,. i i. i 

contributions of I , I , I separately m .big. |5| as a dashed- line, dotted-lme and solid 

curve respectively, we see that each term produces an enhancement near the anti-specular direction. The 

terms I , I are the classic one(see RefsiO) where the fields do not interact with both rough 

(20—20) 

surfaces but produce the peak due to the scattering on the same rough surface: the upper one for I 

and the bottom one for J*' 02 ° 2 '. The contribution J^ 11 11 is the new result of this work, because we 
see that the mechanism of Fig. [I] exists and has the same magnitude as the other terms for the choosen 
parameters. In order to clearly show the displacement of the backscattering enhancement as the angle 
of incidence is varied, we show in Figs. |^ and |?] the numerical results of the perturbation method wj; 
9q = 20°. As expected, we clearly see the peak which is now located at do — —20°. Recent papers 
have also explored the satellite peaks phenomenon which occurs when the wave follows two reverse paths 
but with different guided-mode excitations. In Fig. M, this phenomenon appears for the (TM) to (TM) 

—(02-02) 

polarization due to the term J . Although a similar phenomenon occurs for the other terms it is 

too weak to contribute significantly. 



VII. CONCLUSION 



In conclusion, the results in this paper clearly show that the backscattering enhancement produced by 
a rough slab due to the guided-wave polaritons has several origin. At small order in the perturbative 
development, the peak is produced by a double scattering mechanism. The new results we would like 
to emphasize, is that not only processes where the two scattering events take place on the same rough 
surface are capable to produce the backscattering peak, but it is also due to the fact that the wave can 
be scattered first by one of the rough surface and second by the other surface. To carry the proof, we 
have performed a perturbative development up to order four in the rms-height of the surfaces which has 
been possible using an integral equation called reduced Rayleigh equation, in which only the scattering 
matrix of the whole structure is unknown. 
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APPENDIX A: SCATTERING MATRIX COEFFICIENTS 



Here we collect the expression of the matrices not given inl 



2 l 1 



S ± (P\P ) 



£l - £o 



2a (p )( eo ei) 1 / 2 
After some calculations we obtain : 



^ r - F l-,0+, . . —(00). . . , -=-f1-.0-. i . 



(Al) 
(A2) 



S + (P\P ) 



"(eoei) 1 / 2 



£ i IN \\Po\\ F v(Po) + e a 1 (p)a 1 (p )F v {p Q )p-p ^ K o «i(p) F+(p ) (p x p ) 

( eo ei) 1/2 tf 2 ^(Po)p-Po 



^ (p|Po) = 



^ -eoei /2 ^oai(Po)^y(Po)(P x Po) 

■[oloiPoT 1 , 

(ei - e o) 
ao(Po) (eo ei) 1 / 2 

/ -eoai(p )||p||||p ||F v 7(p ) -e^ 2 K a t (p) ai(p ) ^(Po) (p x P )A 
x -eiai(p)ag(p ) J F^(p )p-p 

Vei /2 if ag(p ) J F 1 +(p )(pxp ) 2 -(e e 1 ) 1 /2 j ^2 Ql ( Po )F H (p )p-p / 

•^(Po)r 1 . 



(A3) 



(A4) 
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FIG. 1: New mechanisms responsible of enhanced backscattering. The incident wave is transmitted in A 
as perfect interface between the media and 1, then the wave is scattered by the second rough surface 
in B and by the first one in C. But the wave can now follow this path the other way round. The phase 
difference between these two waves is zero in the anti-specular direction which produces the peak. 




FIG. 3: Definition of the scattering vector. 
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FIG. 4: The bistatic coefficients for an horizontal (TE) and a vertical (TM) polarized normally incident 
light of wavelength A = 632.8 nm, on a slab with an upper two-dimensional randomly rough surface, 
characterized by the parameters (J\ = 15 nm, li = 100 nm, and a bottom rough surface characterized by 
(i2 = hnm, I2 — lOOnm. The dielectric constants are ei = 2.6896 + i 0.0075 and €2 = — 18.3 + z0.55. 
The thickness of the film is H — 500 nm. The scattered field is observed in the incident plane. For 
each figure are plotted : the total incoherent scattering ^ ncoh (solid curve), the first order given by 
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FIG. 6: Same parameters as in Fig. [|, but with 6q = —20' 
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FIG. 7: Second order contributions to Fig. ^ 



